Phonon and Electronic Non-radiative Decay of Excitons in Carbon Nanotubes 
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We investigate theoretically the rates of non-radiative decay of excited semiconducting nanotubes 
by a variety of decay mechanisms and compare with experimental findings. We find that the multi- 
phonon decay (MPD) of free excitons is too slow to be responsible for the experimentally observed 
lifetimes. However, MPD lifetimes of localized excitons could be 2-3 orders of magnitude shorter. 
We also propose a new decay mechanism that relies on a finite doping of nanotubes and involves 
exciton decay into an optical phonon and an mtraband electron-hole pair. The resulting lifetime is 
in the range of 5 to 100 ps, even for a moderate doping level. 



Semiconducting single-walled carbon nanotubes 
(CNTs) are one-dimensional, direct band-gap materials 
with strongly-bound exciton states, that are actively 
being explored for technological applicationsi. The 
potential usefulness of CNTs as an optical material 
depends, to a large extent, on their luminescence 
efficiency, which is determined by their radiative and 
nonradiative lifetimes. The CNT radiative lifetime has 
been calculated^i^ to be in the range of 1-10 ns. The 
measured lifetime, however, is much shorter, in the 
range of 10 to 100 psii^i^iiiSii^, indicating the presence of 
an efficient non-radiative decay channel(s). The nature 
of the dominant non-radiative decay path has been the 
subject of debateiSiiiii^. Low energy, 'dark' exciton 
states^i^ii^ii^ have been proposed as the cause for the 
small yield, but the small dark-bright exciton splitting 
does not support this conclusions^. It is also known 
that excitons in bulk semiconductors can be localized 
at structural defects, impurities, grain boundaries and 
other heterogeneities^^. Analogous exciton localization 
can occur in CNTs. In fact, the quasi-lD, single atomic 
layer structure of single-walled CNTs makes them quite 
sensitive to potential fluctuations in their environment, 
e.g. trapped charges on an insulating substrate, adsorbed 
species, etc. Recent photovoltaic studies of CNTs have 
directly revealed such potential fluctuationsi^. Even in 
solution, the coating of surfactant-coated CNTs may not 
be homogeneous, leading to potential fluctuations, and 
in otherwise perfect tubes their ends can act as exciton 
traps. Photoluminescence studies of surfactant-coated 
CNTs have indeed provided strong evidence for the 
formation of localized cxciton a^^i^^'^^ . The strong effect 
of the environment on the non-radiative decay is most 
clearly shown by the fact that no fluorescence was 
detected from CNTs on an Si02 substrate^S, but the 
tubes luminescence strongly when suspendedSl. It is 
also known that, although not intentionally doped, 
CNTs on most substrates in air behave as p-type 
semiconductorsii. Suspended CNTs, on the other hand, 
usually show intrinsic behavior—. Furthermore, CNTs 
on insulators show pronounced hystcretic effects due 
to shifts of their Fermi level (field doping) by trapped 
charges^^iSiiS^. Analogous charge-transfer doping effects 



have been suggested for a variety of adsorbed gases on 
CNTs^^. 

It is clear that a variety of factors and environmental 
influences can affect the non-radiative decay of excited 
nanotubes. In an effort to evaluate their contributions, 
we present in this Letter a theoretical study of the ef- 
ficiency of non-radiative decay pathways of nanotubes 
involving purely multi-phonon decay as well as electronic 
decay mechanisms. We examine the decay behavior of 
free-exciton nanotube states, trapped (localized) exci- 
tons, and excitons in doped nanotubes, i.e. nanotubes 
whose Fermi level does not coincide with the neutral- 
ity level. In the latter case, we propose and evaluate 
a novel, very efficient, phonon-assisted indirect exciton 
ionization (PAIEI) mechanism. Our calculations indi- 
cate that the combination of localized exciton MPD and 
the PAIEI mechanism allow us to explain the range of 
available experimental data on the non-radiative lifetime. 
Furthermore, the calculations suggest that, in principle, 
a high emission quantum yield can be achieved in intrin- 
sic, small diameter tubes. 
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FIG. 1: Schematics of the exciton decay mechanisms: (a) 
Multiphonon Decay (MPD) and (b) Phonon- Assisted Indirect 
Exciton Ionization (PAIEI) in p-doped CNTs. 

First, we explore the MPD decay process for free exci- 
tons schematically shown in Fig. [T^. The MPD decay of 
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excited CNTs to the ground state has been proposediS, 
as being plausible for CNTs because of the high energy 
C-C stretching phonon (G-mode, hoj w 0.2eV) and its 
strong electron-phonon coupling2ii2^. 

To calculate the MPD rate we need to evaluate the 
nth order diagram depicted in Fig. [1^?^ . The electron- 
phonon Hamiltonian involves both intraband and inter- 
band (Af^) electron-phonon couplings. It is obtained us- 
ing the Su-Schreiffer-Heeger model as in refill including 
a bond-bending coupling term^Sii^S. The exciton-phonon 
coupling is derived from the electron-phonon coupling 
using the two-particle exciton wavefunctions calculated 
from the Bethe-Salpeter Equation as in refj ^'^i^^ . 
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FIG. 2: (a) MPD rate in four nanotubes as a function of 
exciton energy calculated with e = 3.3 for the four lowest en- 
ergy bands: bright (red), dark (blue), doubly degenerate dark 
(green). The vertical dashed lines show characteristic ener- 
gies of different combinations of the zone-center G mode and 
even numbers of the zone-boundary D modes, (b) averaged 
MPD rate in 118 tubes as a function of tube diameter (red 
circles). The solid red curve is a fit to Eq. ^ reduced by 
n\ with 7 = 0.04 eV. The black curve is a simulation of the 
localized exciton MPD rate. 

The results of the MPD rates for the first four low- 
est energy free excitons are shown in Fig. [2^. When 
the exciton energy matches an integer number of opti- 
cal phonon quanta, the MPD rate shows a strong res- 
onance, see Fig. The Su-Schreiffer-Heeger type cou- 
pling does not break the symmetry of the A and B carbon 
atom sub-lattices in the CNTs and, therefore, does not 
allow coupling among even (dark) and odd (bright) parity 
excitons^. The bond-bending tcrm^ has a much smaller 



phonon coupling, but it facilitates the decay of the odd 
parity states. The bright exciton {L = 0) can decay only 
by emitting phonons whose momenta sum up to zero. 
This is achieved by a combination of even numbers of 
zone-boundary optical phonons (D-modes of about 0.18 
eV) with opposite momenta and arbitrary numbers of 
Raman active G-phonons with zero momentum. 

Considering the limit of a localized single exciton state 
coupled to a single phonon mode, the MPD decay reduces 
to the classical Franck-Condon problem, which is exactly 
solvable in this casei^: 
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where E is exciton energy, Sj and S are the interband and 
intraband Huang-Rhys factors^, and n is the number 
of emitted phonons. The phonon density of states and 
phonon lifetime are modeled by the broadening param- 
eter 7. The interband exciton-phonon coupling changes 
the number of vibrational quanta by one, while the in- 
traband coupling determines the overlap between vibra- 
tionally excited state with n — 1 phonons and vibrational 
ground state. 

The main difference between the calculated MPD rate 
and the rate estimated according to Eq. ^ is the sub- 
stantially smaller magnitude of the former. In CNTs, 
in the emission process of a phonon with momentum 
q' by an exciton with momentum q, the phase of the 
exciton-phonon matrix depends on both q' and q. We 
find that cross-terms corresponding to the phonon emis- 
sions in different orders, leading to the same final state 
with n phonons, will sum to nearly zero, reducing the 
total MPD rate by nearly a factor of nl^. 

We have calculated the MPD rate for all chirality semi- 
conducting tubes in the diameter range from 0.8 nm to 
2.0 nm by averaging the energy dependent decay rates 
in Fig. [5^ over the exciton distribution. The results are 
shown in Fig. ^jp. We find that only a few large diame- 
ter tubes, where three or less optical phonons are needed 
to annihilate the exciton, have lifetimes in the range of 
100 ps to 1 ns, whereas most of the other tubes show a 
much longer non-radiative lifetime. We can also fit the 
calculated rates in Fig. [2j3 to the localized limit expres- 
sion Eq. ID) reduced by n\. The diameter scaling of the 
intraband Huang-Rhys factor is similar to that of the op- 
tical phonon sideband intensit}»21 S oc 1/d, where d is 
a tube diameter. The interband Huang-Rhys factor for 
the free e-h pair scales inversely with the length L of the 
tube Si oc \M^\^/N = A\M^\'^ /{dL), where we calculate 
w 0.1 eV^ and N is the number of two atoms prim- 
itive unit-cells of area A w 5.39 in the super-cell. For 
a bound exciton, we find that the interband coupling is 
enhanced by a factor of L/tq, where the exciton radius 
rg depends on the diameter and effective dielectric con- 
stant e that accounts for screening by both the nanotube 
itself and the immediate environmcnl^. The choice of 
e = 3.3 reproduces the experimentally measured exciton 
binding energy in CNTs^, and we find rp = 1.95(i and 
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Si « 0.02 /(f. Therefore, we are left with just one fit pa- 
rameter S and the best agreement is found for S = O.lA/d 
in Fig. Eb. 

To model localized exciton MPD decay, we assume that 
the intcrband Huang-Rhys factor Sj is the same as that 
calculated for the free excitons. We then calculate the 
intraband Huang-Rhys factor enhancement^ by forcing 
the four lowest energy excitonic bands to have a flat dis- 
persion, while keeping the strength of the exciton-phonon 
coupling constant. Depending on the tube diameter, we 
find an effective S enhancement by 1.4 to 1.8. Using 
Eq. (HD with S = 0.21/d, we obtain the MPD rate of 
localized excitons shown in Fig. [2b. We find that the in- 
creased effective exciton-phonon coupling, due to the re- 
duced exciton dispersion, and the coherent contributions 
of matrix elements bring the MPD rates for the localized 
states close to the experimentally measured non-radiative 
decay values, even for small diameter CNTs. 

We now examine electronic non-radiative decay mech- 
anisms. A purely electronic (Auger) decay mechanism^ 
can be viewed as the reverse of an impact excitation 
process^. However, the angular momentum conserva- 
tion law sets strong restrictions for the Auger decay, such 
that the En bright and dark excitons with zero angular 
momentum cannot decay through such a process^. The 
higher energy dark exciton with finite angular momen- 
tum can decay, but only at a high doping level with the 
Fermi level lying above the bottom of the second band^^. 
We also note that at high excitation densities, exciton- 
exciton annihilation dominates the decay^i^^. However, 
even at low excitation densities the experimentally mea- 
sured decay is fast. 

We show here that in the presence of the free car- 
riers in CNTs, an exciton can decay fast by creating 
a phonon and an intraband electron-hole (e-h) pair, as 
shown in Fig. [T|3. This phonon-assistcd indirect exci- 
ton ionization (PAIEI) process involves Coulomb and 
exciton-phonon interactions which are both very strong 
in CNTsHi^iS^i^^iSl. The interband exciton-phonon in- 
teraction allows the exciton to cross to the electronic 
ground state. However, emission of a single phonon is 
not sufficient to conserve energy. The excess energy is 
accommodated by producing an e-h pair in the valence 
band (for a p-doped CNT). The initial excitonic state 
is described by a two-particle wavefunction of exciton in 
state p with momentum j^I/P >= ^fcg'^I+g^fcl^'S' >• 
Here, I^^^P gives the contribution of the final hole state 
k with energy Ek to the superposition of states describing 
the initial exciton state. The PAIEI decay rate is given by 
Fermi's golden rule as shown by the diagram in Fig.[T]3i2S 
and is controlled by the overlap in reciprocal space of the 
excitonic and final hole state wavefunctions I^^^P, and 
by the Fermi level which determines the available final 
states for scattering. 

The computed energy-dependent PAIEI rates are 
shown in Fig. [5^ for the four lowest energy excitons in 
several nanotubcs for a fixed doping level of p = 0.21 
e/nm. Unlike in the MPD case, which shows resonances 
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FIG. 3: PAIEI decay rate: (a) as a function of exciton energy 
for the four lowest energy bands: bright (red), dark (blue), 
doubly degenerate dark (green); (b) averaged PAIEI decay 
rate as a function of free carrier density in the same tubes: 
red - (13, 0), blue - (16, 0), green - (19,0), and cyan - (25, 0) 
using £ — 3.3; (c) averaged PAIEI decay rate in 118 tubes at 
constant carrier density p — 0.21 e/nm, using e — 2.5 (blue 
circles), e = 3.3 (red circles), e — 4.0 (green circles) along 
with the simplified expression Eq. ((2)1 (solid lines). 



for an integer number of optical phonons, the PAIEI rate 
shows a much weaker energy dependence. We find that 
the PAIEI rate at the bottom of the each band in all 
chirality tubes is largest for the lowest energy dark exci- 
ton, while the finite angular momentum, higher energy, 
dark exciton has the smallest PAIEI rate (Fig. [3^). As 
the doping increases, the PAIEI rate increases propor- 
tionally (Fig.H)). 

In the low doping limit, which is the main interest 
here, the EIPAI rate can be found analytically. The num- 
ber of final electron states is proportional to the Fermi 
wavevector kp, which defines the free carrier density p. 
At zero temperature, kp = irp/A. This explains the lin- 
ear dependence of the PAIEI decay rate on p, which 
holds up to a 0.4 e/nm doping level. The number of 
final hole states is given by the density of states, which 
we calculate using hyperbolic band dispersion^, at the 
hole energy Sh ~ 3A — hui, where A « 0.42 eV nm/d 
is half the bandgap. The tail of the exciton wavefunc- 
tion in momentum space is the key factor determining 
the PAIEI rate. In the case of an uncorrelated (free) e- 
h pair, the two particle wavefunction is a delta function 
and the PAIEI decay cannot take place. We find that 
the reciprocal space wavefunction of the lowest energy 
exciton is best approximated by the following functional 
form |A^=LoP oc (1 -h fcV^)-^-^.^. The excitonic radius 
To is a fit parameter here and scales with tube diameter 
d and effective dielectric constant e as ro « 0.824de° ''^. 
The analytical expression of the PAIEI decay can then 
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be written as: 

1 _2TT ASih^u'^ khrp p 
TpAiEi h 2TrCd (^i + k2rlfe s,,' 

where kh = ^Je\ - lS?/(hVF), Vp « 10*' m/s is the Fermi 
velocity of graphene, C « 0.65 is the wavefunction nor- 
mahzation constant, Si is an effective interband Huang- 
Rhys factor. Eq. ^ fits remarkably well the PAIEI 
rate at room temperature and fixed carrier concentra- 
tion for all chirality semiconducting tubes with diameter 
range between 0.8 nm and 2.0 nm and for a variety of e 
with just one adjustable parameter 5*1 = 2.5 (Fig. [Sh). 
This value of 5*1 suggests an electron-phonon coupling 
|M^P = Sih^Lo"^ = 0.1 eV^, which is in good agreement 
with the calculated value here and measured in^^. Fi- 
nally, we also note that the PAIEI decay rate decreases 
with increasing effective dielectric constant due to the 
increased exciton radius. 



In summary, we have calculated the non-radiative life- 
times of excited semiconducting CNTs via multiphonon 
decay and by a new, phonon-assisted electronic decay 
channel, function of the CNT diameter. We find 

that the MPD rate for free excitons is too slow, espe- 
cially for small diameter CNTs, to account for the exper- 
imentally measured decay rates. On the other hand, we 
predict that the MPD decay rate of localized excitons is 
increased by a few orders of magnitude, due to a coher- 
ent contribution of exciton-phonon matrix elements and 
a stronger effective exciton-phonon coupling. Thus, lo- 
calized exciton MPD rates can account for a wide range 
of available experimental data. Exciton localization is 
usually associated with charge transfer and doping. In 
these doped tubes, we find that the decay channel in- 
volving phonon-assisted indirect exciton ionization can 
dominate both the free and localized exciton decay, even 
at moderately small doping levels. 
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APPENDIX A: AUGER DECAY 



The Auger decay of the exciton into electron-hole pair involves Coulomb interaction, which is very strong in carbon 
nanotubeaii^. In fact, Auger decay^ is a reverse process of the impact excitation, which is much faster in CNTs than 
in bulk semiconductors^. However, for the Auger decay to take place both energy and momenta have to be conserved: 

where E is the exciton energy and e is the electron and hole energies. The lower index denotes the continuous 
momentum along the tube axis, and the upper index the discrete angular momentum, which can have values I ~ 
±1,±2,±4 etc. The angular momentum conservation law sets strong restrictions for the Auger decay. For example, 
the bright En exciton with zero angular momentum and q = can not decay by Auger. Indeed, the right hand side 
of Eq. I|A1[) is identically zero for L = and q = 0. A similar argument applies for the lowest energy dark exciton with 
AL = 0. The higher energy dark excitons^ with finite angular momentum AL = 2 can decay by Auger, but in order 
to conserve the angular momentum, a hole has to be created in the second band = 2 and an electron in the third 
band 1^ = 4. This requires a high doping level, so that the Fermi level lies above the bottom of the second band. The 
lowest energy exciton which can decay by Auger at low doping level is E12 with AL = 3, which can be photoexcited 
with light polarized perpendicular to the tube axis. Indeed, both momentum and energy can be conserved when E12 
exciton decays into hole in the first band Z/i = 1 and electron in the third band = 4 as shown in Fig. [H 





FIG. 4: Schematics of the Auger decay, which can take place at low doping level for a E12 exciton with angular momentum 
AL = 3. The Fermi level Ef is shown by the dashed line. The angular momentum conservation prevents the optically bright 
exciton E\\ {AL = 0) to decay by this mechanism. The higher energy, dark excitons with AL — 2 can decay by the Auger 
mechanism, but only at a high doping level with Ef lying above the bottom of the second band. 



APPENDIX B: THE MPD RATE CALCULATIONS 



In this section we derive the multiphonon decay (MPD) rate equations using perturbation theory expansion. The 
transition rate from an initial state i to a final state / is given as: 

W.f = ^S{Ef-E,)\{f\mf 
/f\T\-\ /.lT.l■^ , , {f\V\k') {k'\V\k) {k\V\t) , 

where V is the exciton-phonon potential and 7 is the level broadening. 

In our case, the initial state exciton two-particle wavefunction with momentum q in state S is given by Eq. (|F2[) 
and the final state |GS, 0;^^^^ , Wq^^j , w^g^g , ...) is the electronic ground state |GS) with a number of created phonons 
with all possible momenta q and band indexes fi during the decay process. The number of created phonons in the 
final state determines the order of the perturbation. 
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The total MPD rate of an cxciton {q, S) is given by the sum of the first, second etc. orders of perturbation: 



qS 



qS 



(B2) 



i=l 



The lowest order non-zero contribution determines the MPD rate. The effective decay rate is averaged over the exciton 
distribution. In the case when the thermalization is fast compared to the non-radiative lifetime, it is given by: 



Wn 



exp 



kf^T 



(B3) 



where index v = [q, a, S) includes summation over exciton wavevector q, spin index a (labeling singlet and triplet 
states), and exciton state S (labeling the bound and continuum states for a given q and a). The partition function is 
Z{T) = exp{—E^/kBT), summed over the same range of v. 

Below we give explicit expressions for the first, second, and third order decay rates: 

First order (Fermi's Golden rule): 
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Second order: 
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Third order: 
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(B6) 



where Pqiq2q3 stands for the sum over six permutations of indices qi, q2, and q^. The interband B^q,^ and intraband 
D^^ exciton phonon couplings are given in appendix |F1 

Note, that in the absence of exciton and phonon dispersions, Eq. (|B4[) - (|B6[) reduces to the MPD rate in the 
localized limit given by Eq. (2) of the main text, except for the missing factor e^"^ « 1, which can not be obtained by 
perturbation theory. 

We find that the cross terms in Eq. (|B6p . corresponding to the different orders of emitted phonons, sum nearly to 
zero. Therefore, neglecting the interference between these terms is a good approximation to the total MPD rate: 
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The evaluation time of Eq. (|B7p scales as N" ^ which becomes computationally prohibitive for n > 4. Alternatively, 
this expression is equivalent to the following set of equations: 



FiieqSS') 
FNieqSS') 



qifj.1 SiS{ ' 



FN{E^,q,S,S) 
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(B8) 



which are used to calculate the MPD rates in Fig. 2 and 3 of the article. 

In Fig. [5] we calculate the MPD rate using the exact expression, as in Eq. (jB4p - (|B6p and neglecting interference 
terms, as in Eq. (|B8[) . The coherent terms contribute about 1.5 times to the MPD rate for n = 3 and about 2.5 times 
for n = 4, which are significantly lower than in the fully coherent case which give a factor of 6 larger rate for n — 3 
and a factor of 24 for n = 4. 

In the test calculations we replaced the exciton-phonon matrix elements in Eq. (|B6[) by their absolute values (forcing 
their phases to be identically zero) and found that the MPD rate for order n = 3 increases by a factor of 4, which is 
much closer to the fully coherent result of 3! = 6 times larger rate. 
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FIG. 5: (a) MPD decay rate of free excitons as a function of tube diameter including coherent contributions (green circles) and 
neglecting them (red circles), (b) The ratio of the two calculations (black circles) with a two dashed lines corresponding to 1.5 
times and 2.2 times rate enhancement due to the coherent contributions for n = 3 and n = 4 correspondingly. 



APPENDIX C: PAIEI DECAY RATE AND EXCITON WAVEFUNCTION 



In this section, we give the analytical form of the exciton wavefunction in the reciprocal space, which determines 
the rate of the Phonon- Assisted Indirect Exciton Ionization (PAIEI) . Exciton in state p with momentum q and energy 
E'P can PAIEI decay with the rate given by the Fermi's Golden rule: 

^ ^ H \^^k^+qM-k^-qJ^\^U^^ {El - hujq+k,-k, " £fci + £fc. ) (1 + nq+k^-k,,^.) fkr (1 - fkj (CI) 

TPAIEI ^"^k.k^f. 

where the Bose-Einstein factor Ug gives the phonon population, while the Fermi-Dirac distribution fi; defines the 
probabilities that the hole state with energy and momentum (ej.j,/ci) and the electron state {£k2T^2) in the valence 
band are available for scattering. 
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The exciton wavefunction A^^ gives the contribution of the final hole state to the superposition of states describing 
the initial exciton state. We find empirically that it is best fitted by the following equation: 



N\Ai\' = -— ^^^^^ (C2) 



with just one adjustable parameter the exciton radius rg. Here, the normalization constant C = dx/ (1 + x^)^'^ « 
0.65 requires l^f^P = Ij and A = 0.0539 nm^ is the area of the two-atom primitive unit cell, d is the tube 
diameter, and N is the number of primitive unit cells in the supercell. The results of the fits are excellent, as shown 
in Fig. [6] for few selected tubes. We find an empirical scaling of the exciton radius with tube diameter and effective 
dielectric constant e, that accounts for screening by both the nanotube itself and the immediate environment, to be: 

ro = 0.824 d e°-^2 (C3) 

which is similar to the scaling found in reff^ from the real space wavefunction analysis. 




FIG. 6: The lowest energy excitonic wavefunctions in reciprocal space in three tubes (13,0) - red circles, (19,0) - green circles, 
and (25,0) - red circles calculated for e — 3.3. The black solid curves (shown behind the circles) are best fits to Eq. (|C2|I with 
ro=2.05 nm, 2.71 nm, and 3.36 nm, correspondingly. 



APPENDIX D: SSH ELECTRON-PHONON MATRIX ELEMENTS 



In this section we give explicit form of the Su-Schreiffer-Heegcr (SSH) Hamiltonian^i^, in reciprocal space. The real 
space SSH Hamiltonian reads: 



iS 



XI cxp {ikR^+s 



^H^g exp (zgT?,) 



N 

9 



n^s = rf+s-rf, (Dl) 

where i?j is the position of the unit cell i; rf is the position of atom A in unit cell i; N is the number of two atom 
unit cells in the supercell; g = 5.3 eV/A is the electron-phonon coupling^. 
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The SSH Hamiltonian in reciprocal space is: 



kq6 

X {ctl.C^ exp (tk (i?,- - i?o) ) + Cl^l.C^ exp (-t (k + q) - ^o) ) ) , (D2) 

where is the rotation operator acting on displacement vector in the i = unit cell and rotating it to i = (5 unit 
cell. In the tight-binding basis: 



Ik 



^exp(-^fc(i?^-i^o)) (D3) 



5 



the electron-phonon interaction has the form: 



kqfj, 



where 



Ml,^ = M^^ = L exp Uq(Rs -Ro))-A 



2 Y 2Muj^g ^ - 

X ( exp (ik (Rs -Ro))+ cxp {-i{k + q] [Rs- Ro 



7fc| ^ ^ I7fe+g 



X 



exp (ik - i?o) ) - exp (-Z (fc + - Ro))j (D5) 



The phonon displacements are written in the second quantized form: 



^9 



where Ag^ and are amplitudes of the vibrations with frequency oj^, and M is carbon mass. 
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APPENDIX E: BOND-BENDING ELECTRON-PHONON MATRIX ELEMENTS 



111 this section we give explicit form of the bond-bending Hamiltonian in reciprocal space. The bond-bending 
electron-plionon Hamiltonian in real space is: 

i6i ^62 iSi ^62 



. , /r -c N COS {(j))nAi6i - nAiS2 , cos {(j))nAi62 - nAiSi 

A(j)A {01162) = u,^s + — — u,^s^ 

VAiSi sm (0) rAi52 sm (0) 

, (rAiSi - rAiS2 cos (</>)) n^ia^ + (r^i^a " r MSi cos {(f))) fiAi52 -^a 



riAiS = 



rAi5irAiS2 sin(0) 



\rAis\ 



cos((/)) = nAiSi'riAt52 (El) 

where A0 is the C-C-C bond angle change, Ri is the position of the unit cell i; rf is the position of atom A in unit 
cell i; N is the number of two atom unit cells in the superccU; gi, ~ 0.9 eV/rad is the clectron-phonon coupling^. 

The bond-bending Hamiltonian in reciprocal space reads: 



nfi;i = 7= E Di{qSiS2)C;^l^C^ exp [ik (^Rs2 -Ro)-i{k + q) (i?^, - R 

V kq5i^52 

+ 7^ E ^2 {qSiS2)C^lgC^ exp (ik -Ro)-i[k + (^Rs, - R 

V ^ kqSi^S2 

COS {(f))nASi ~ nAS2 frf. ^bA --/n u \ I cos{(j))nAS2-nASi (rt. .^/^ 5 

TASi sm (0) V / V / rAS2 sm (</>) V / V 

^ - '•^^2 cos ((/>)) n^ji -I- (r^^a - TASi cos {(/))) nAS2 ^A 

?'A5i''A52 sin(0) 

r> f X X \ cos (0)nBg, - ng^, /'.f, --/^ g \ , cos {(f))nBS2 ~ nss^ frf, -^a\ -^fff „ 

D2(q5id2) = T^jU cxpzg - i?o H ^-t-t T^, it exp i?^^ - i?o 

rB5i sm (0) V / V / 7-552 sm (0) V / V 

, (rs^i - rBS2 cos (0)) ns^i + (rsia - cos (0)) nBS2 /T^o^ 

''B5i?'SA-2 sin(0) ^ 

where T5 is the rotation operator, which rotates carbon atom in the original unit cell to the carbon atom in the 
neighboring unit cell 6. In the tight-binding basis Eq. (jD3[) the bond-bending coupling is: 



kqfi 

where 



Ml,^ = -^4',;. = F2 {kqy) - Fi (fcg/i) (E4) 



11 



where 



Fi{kqn) = ^ ^ exp (ik (Rs^ - Ra) - i (k + q] (Rs^ - Rq 



V dl/()2 

X ^-7-: Ts^ Bg^ exp iq Rs^ - Ro 

rASi sm (0) V y V 

cos {(j))riA52 - riASi frf. g \ .^/^ g 
H T^a-B^M exp iq Rs., - Ro 

, {taSi - 1'A52 cos ((/>)) nASt + (T'Aga ~ ^^A^i COS {(/))) UaS^ 7* 

rASirAS2 sm (0) 



V ' 5i^S2 
COS {4))nBSi - ?IB,5 



TsiAqfj,j cxpiqyRs-^ - Ro 
fs^Aq^^ cxpiq(^Rs2 - Ro 



VBSi sin (0) 

cos {(j))nBS2 ~ ^-B.5i 

rBS2 sin (0) 

(^•Bai - rBS2 cos (0)) ngg^ + (rg^^ - rg^^ cos ((/))) nBS2 
TBSirBSi sin((/)) 



S.m] (E5) 



APPENDIX F: EXCITON-PHONON MATRIX ELEMENTS 

In this section we derive the exciton-phonon Hamihonian starting with the electron-phonon Hamihonian: 

Hel-ph =^Y1 {^'^lqt^^l+q'"k ~ Mlq^^cl^gCk + Mfc^^cj^+^ffe - Af^^^ V^^^ Cfe ) (fl,^ + oL^^) , (Fl) 



where M^^^ (x g is momentum dependent electron-phonon coupling (see appendixes ID] and IE)) : is a phonon 

creation operator with wavevector —q and phonon band index fi = 1...6; and N is the number of primitive unit cells, 
each containing two carbons; cl_^^ (vk) creation (annihilation) of an electron in the conduction (valence) band acting 

on the ground state |GS) = Ylk "^t |vac). The indices k and q label both the continuous ID wavevector along the tube 
axis and the band index. 

To derive the exciton-phonon matrix elements from Eq. (jFip we need to know the exciton wavefunction given by 
the solution of the Bethe-Salpeter Equation (BSE) for the two-particle exciton wavefunctioni*^: 

fe 

Here A^^ is the eigenvector of the S's state of BSE solution. Using the orthogonality relation of the BSE solution 
(Es^'ig^i*q = ^kk'), we obtain: 

cl+,^;fc = ^<|*f). (F3) 
s 

Using Eq. (jFHIF3[) . we derive the exciton-phonon Hamiltonian: 



g' 



^ S'q'ii 

(v/"9M + 1 + 1) + xA^h-gp - 1)) |GS) 

fe 
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where Ug^ is the Bose-Einstein phonon occupation number. Here we neglect the couphng of the one exciton subspace 
to the subspace with two excitons. Matrix elements B^^i^ are the intraband and D^^ are the interband exciton-phonon 
couplings. 
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